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Recently a number of integrals and Euler sums of interest in Physics and Quan- 
tum Chemistry have been shown to be interrelated [T] . When calculating radia- 
tive corrections in Quantum Field Theory, multi-dimensional Feynman integrals 
which lead to these sums are encountered. For typical applications to Quan- 
tum Field Theory, see Cheng and Wu [2]. A common theme in these works 
has been the presence of Neilsen's generalized polylogarithm functions [3] . In 
' this article we continue these investigations and provide a number of new iden- 

tities involving integrals of logarithms, polylogarithms, Euler sums and Neilsen 
functions. 



2 Polylogarithm Integrals 



Integrals containing the product of two polylogarithmic functions Li r (t) [4] of 
orders p and q and having the forms (p,q > 1) 



, tM ./' 5M* . w 
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J_(p,9)= f 1 Li P (-t)ZA q (-t)dt =I _ (qt 



*,P), (2) 

J ± (p,g)= f L ^ Ll ^ d \ ( 3 ) 

arise in a number of physical applications [3] and are of interest in themselves. 
The integral I + (p, q) has been studied by Frietas [6] . We note the p, q symmetry 
displayed in Eqs. (1) and (2). Integration by parts for any of the integrals above 
allows one to obtain the following partial difference equation 

I(p,q)+I(p+l,q-l) = R{p+l,q). (4) 

A more symmetric form of this difference equation being 

I(p,q-l)+I(p-l,q) = R(p,q), (5) 

where I(p,q) stands for any of the integrals in Eqs. (1-3), and R(p,q) stands 
for one of the three corresponding quantities 

R + (p,q)=Li p (l)Li q (l) = C(p)((q), 

R-(p, q) = Li p {-1) Li q (-1) = (2 1 -? - l)C(p) (2 1 -* - l)C(g), 

R±(p, q) = Li p (l)Li q (-l) = C(p) - 1) C(?), 

and C(z) is the Riemann zeta function. We note that in the special case where 
z approaches 1 that 

hm(2 1 - z -l)C(z) = -ln(2). 

z— >1 

A set of solutions for the difference equation Eq. (5) is given by 

I(p,q)-J2(- 1 ) k R(P + k + l,q-k) . (6) 



I(j) + n,q- n) = (-1)" 



fc=0 

Proof : If p is replaced by p + k and q is replaced by q — k in Eq. (4) we get 

I{p + k,q - k) + I(p + k + l,q - k - 1) = R(p + k + 1, q - k). 
Multiplying this by (— l) k and summing k over n terms produces 

n — 1 n — 1 n — 1 

^(-l) fc I{p+k, q~k)+J2(-l) k I(p+k+l, q-k-l) = J](-l) fc R(p+k+l, q-k). 
k=0 k=Q k=0 

Writing out the last term in the second sum and the first term in the first sum 
on the lhs of this equation produces 

n-l 

(-l)"- 1 /^ + n,q- n) + I(p, q) = ]T(-l) fc R(p + k + l,q - k), 

the remaining sums having cancelled. Rearrangement of this equation gives the 
desired result. 



2 



3 Special cases of the I(p, q) integrals 



3.0.1 The I + (p,q) and I_(p,q) integrals 

The odd case q — p = 2n — 1 In the case of I+(p, q), and I-(p, q) we 
note that when the integrals' parameters q and p are set equal and substituted 
in Eq. (5) that I(p,p — 1) = \R(p,p), or more usefully 

I(p + n,p + n- 1) = \R{p + n,p + n), (7) 

i.e. the I integrals can be written in closed form. This follows directly from the 
p, q symmetry of these integrals and Eq. (5). Consequently, the family of all 
such integrals where q and p differ by an odd integer can likewise be expressed 
in closed form. To see this we set q = p + 2n — 1 in Eq. (6) and solve for 
I(p,p + 2n — 1) having used Eq. (7) and noting that the last term in the sum 
also contains the quantity R(p + n,p + n). Wc find 

n-2 

I(p,p+2n-l) = ±(-l)" +1 i?(p + n,p + n) + ^(-l) fe R(p+k + l,p + 2n-l- k), 

k=0 

(8) 

that is, these integrals are given by a finite sum of zeta functions. Using the p, q 
values associated with the series defined by p + q = 5 we have for either I + (p, q) 
or I_(p,q), the examples 

7(l,4) = -ii?(3,3) + i?(2,4), 

7(2,3) = ii?(3,3). 

The even case q — p = 2n In the case where q and p are set equal in 
Eq. (4), the integrals I + (p,q) and I-(p,q) have the property 

I(p + l,p - 1) = R(p + l,p) - I(p,p). 

All of the integrals in these families can be related to the ones which contain 
the square of a single polylogarithmic function. Setting q = p + 2n in Eq. (6) 
we get 

n-l 

I(p,p + 2n) = (-1)™ I(p + n,p + n) + Y / (-l) k R(p + k + l 7 p + 2n - k). (9) 

fe=0 

By way of example, for the series associated with p + g = 6we have for I + (p, q) 
or 7_ (p, q) 

1(1, 5) = 1(3, 3) + R(2, 5) - R(3, 4), 
J(2,4) = -J(3,3) + i2(3,4). 
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3.0.2 The I±(p,q) integral 

The odd case No simple relationship such as Eq. (7) exists for 
I±{PiP — l)i instead there is the slightly more complicated form obtained from 
Eq. (5) where p and q have been set equal and p has been replaced by p + n. We 
get 

I±(p + n,p + n — 1) + I±(p + n — l,p + n) = R±(p + n,p + n). 

Setting q = p + 2n — 1 and substituting I±(p + n,p + n — 1) into Eq. (6) we 
get upon writing out the last term in the sum 

J±(p,p+2n-l) = {-l) n+1 I±{p+n-l,p+n)+^{-l) k R±(p+k+l,p+2n-l-k). 

Here we see that knowledge of the integrals I±(p + n — l,p + n) is required to 
compute I±(p,p + 2n— 1). 

The even case Relations similar to Eq. (9) exist for the I± (p, q) 
integrals. If q and p are set equal in Eq. (4) 

I±{p+l,p- 1) =R±{p + l,p) -I±(p,p). 

The integrals I± (p + 1 , p — 1 ) are seen to be related to I± (p, p) which contains 
two polylogarithmic functions of order p but with arguments differing in sign. 
Setting q = p + 2n in Eq. (6) we get a result which has exactly the same form 
as that for I + (p,q), or I-(p,q) as shown in Eq. (9) i.e. 



n-l 

I±(p,p + 2n) = (-1)™ I ± {p + n,p + n) + Y,{-^) k R±{P + k + l,P + 2n-k). 

k=0 



4 Polylogarithm / integrals with low order 



In cases with low order, some special integrals arise i.e. for q = 0, 1 we have 

--'-<*■». 

Jo l — t 
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Integrals of the sort shown above arise in the calculation of hadronic heavy quark 
production [7] ■ These are examples of multiple polylogarithms Li mi ,..., mn (xi, x n ) 
as defined by 



Lirn\T •'•^m n (^1) **•) *&7i) ^ ^ 



>m„ , m 2 jltti 
k n >...k 2 >k 1 >0 2 1 



Expressing the integrals above in these terms we have 
nl Li p (±t) 



I i^) dt = _ £ip)i(Tl _ 1)j 
Jo 1 + 1 

f 1 Li p {t) - Lz p (l) 

y ^ o* = £i p ,i(i,i) - C(p + i)) 

* ^(-^-^C- 1 ) ^ = Llpl{h _ 1} + (1 _ 2 -p )c(p + 1} . 



5 Infinite series representations for the / inte- 
grals 

As noted above, for an exact calculation of the I + (p,q), I^(p,q) and I±{p-,q) 
integrals, closed form expressions for the integrals I(r, r) are needed for either 
even or odd values of r. In the case of the quantities I±{p, q), integrals of the 
form I±(r,r + 1) are also required. As will be seen below, it is possible to 
express these integrals as infinite series whose sums can be found in closed form 
in at least some cases. 



We will begin with i_ (p, q) . Expanding one of the polylogarithm functions 
say Li q (—t) in a power series within its integrand , we get 

/ -(P'9) = E i W- / Li p {-t)t k ~ l dt. 

k=l ^° 

The moment integrals appearing in this sum can be obtained by repeated inte- 
gration by parts and we have the general result (with empty sum for p = 1) 

£ Li P (~t) t"- 1 dt = t^ mk+1 y m2+ i)]+f2 t^(2 1 ^-l) C(M). 

where ip(z) is the Psi (Digamma) function. Boyadzhiev has computed the mo- 
ments of Li p (t)^\. The integral I-(p, q) then becomes 

J_(p, q) = (-If f^-l)"^-* - 1) C(M) {V l - p - q - 1) C(P + Q + 1 - M) 
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fe=i 

The infinite sums in this expression can be represented in a more convenient 
form by writing out the tp(k/2 + 1) portion of the latter sum in terms of its even 
and odd values of k. We get 

J_ (p, q) = {-iy 2 [ln(2) (2-"-« - 1) C(p + <?) + (!- 2-p-^ 1 ) C(p + g + 1)] 

+ (-if f^-i)^ 1 -" - 1) C(m) (2"- p -* - 1) C(P + « + 1 - m) 



^ (~l) k mk + 12+7] [j^fc + 1) + 7 ] [VKfc + 1/2) - V>(l/2)] 

fcp+g (2/c)P+9 ^ (2fc+l)P+9 

.fe=l fe=l v ' fe=0 v 7 



where 7 is Euler's constant. In a similar way the integrals I+(p, q) and I±(p, q) 
are given by 

V 00 

I + (p,g) = (-l)* ^(-irCMCfplg + l-riK-ir^-^ + ll + T], 

p=2 fe=l 

/±(p,«) = (-i) p £(-1)" cm (2^- p - 9 -i) c(p+g+i-/i)+(-i) p+1 E ^Mfc+i)+7]. 

M=2 fc=l 

Although not manifestly obvious, the p, q symmetry is preserved in the infinite 
sum relations for I + (p,q) and I_(p,q). 



6 Closed forms for the infinite series 



In the expressions for the / integrals given above, four infinite series containing 
Psi functions occur. We will show below that it is possible to give closed form 
expressions for two of them whereas in the remaining two, closed forms can be 
given only for even values of p + q. 

We begin with a discussion of the Euler sum S+(r) of order r i.e. 

00 

k=i k 

A closed form expression for this sum where r > 2 had originally been given by 
Euler i.e. 

r-2 

S + (r) = ±(r + 2)C(r + l)-±£C(M+l)C(r-M)- 
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In the case of the alternating Euler sum S- (r) i.e. 

s-M = E4^ (fc + 1)+7] ' ( 10 ) 
fe=i 

the situation is more complicated. Rewriting the series in Eq. (10) as sums 
corresponding to the even and odd values of k respectively, we find 

~ ty,(2fc + l) + 7 ] f [^(2fc + 2)+ 7 ] 

(2jfc) r 2. (2fc + ir • 

Each of the Psi functions in these sums can be replaced by functions of lower 
argument using the Psi function's duplication formulas. We get 



fc=0 v ' k=0 y ' 

This expression can be simplified with the help of four known sums which we 
call the Jordan sums Ji(r), 3i(r) ; the Milgram sum A4(r), and the sum C(r) 
where 

M^f Wt X/ fll , (Ha) 



fe=i 



(2/c) r 
[V(* + l)+7] 



w 2 ^ {2k + l) r v ' 

C W = 2L — — = — (Hd) 

The sum £_ (r) can then be rewritten as 

S-(r) = J 2 (r) - Ji(r)+C(r) - M(r) - (1 - 2"'- 1 ) C(r + 1). 

Closed form expressions for the sums J\ (r) and 3%{f) have been given by Jordan 
[10] only for even values of r i.e. 

Jx{2n) = -i(l - 2- 2 "- 1 ) ({2n + 1) + ln(2)(l - 2~ 2 ") £(2n) 

n-1 

- 2- 2 "" 1 ^(2^ - 1) C(2 M ) C(2n + 1 - 2/.), (n > 1), 
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and 

n-l 

J 2 {2n) = |(1 - 2- 2 "- 1 ) C(2n + 1) - £ (2~ a " - 2- 2 "- 1 ) C(2/x) C(2n + 1 - 2 M ). 

The sum .M(r) has been given by Milgram [9] for r > 2 in closed form (with 
empty sum for r < 2) as 

M(r) = ±r (1 - 2-^ 1 ) C(r + 1) - hi(2) (1 - 2" r ) C(r) 

- ^ X> + W 2 "* 2 - x ) <^ + 2 )( 2_M_1 - 2 ~ r ) c(r - 1 - M)- 

This Milgram sum can be simplified considerably because of the reoccurrence 
of the zeta products within the sum. We get 

M(r) = §r (1 - 2- r - 1 ) C(r + 1) - ln(2) (1 - 2" r ) C(r) 

(r-4)/2 

-| J] (2^ +2 -l)CG« + 2)(2^ l - 1 -2- r )C(r-l-M), for even r 

M=0 



and 

X(r) = y (1 - 2— 1 ) C(r + 1) - ln(2) (1 - 2^) C(r) - ± 



(1 - 2"( r 2 1 )v ii— -■ 



2 



(r-5)/2 

-3 £ (2^ I+2 -l)C(/i + 2)(2^- 1 -2- r )C(r-l-M), for odd r 

The closed form expression for the sum C(r), has already been given above as a 
multiple of S+(r). 

Finally, the values of the / integrals can be written in terms of these sums 

as 

j+(p, q) = (-irj2(-ir cm c(p + 1 + 1 - M) + (-l) P+1 5 + (p + 5), 

/±(p,g) = (-ir^(-irc(M) (2^-"-i)cb+g+i-M)+(-i) p+1 5-b+ 9 ), 

J_(p, g) = (-l) p 2 [ln(2) (2^"« - 1) C(p + g) + (1 - 2-"-*- 1 ) £(p + g + 1)] 
+ (-If E(-ir(2 1 ^ - 1) CO) (2^ p - ? - 1) C(P + « + 1 - M) 

M =2 

+ (-l) p (p + g) - 2 C(p + g) + 2 (p + g)] . 
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To summarize the situation, we have found that closed form expressions for 
I+(p, q) and I-{p, q) when q = p + 2n — 1 are given by Eq. (8). In the cases 
where q = p + 2n the sums 5+(2p + 2n) and S-(2p + 2n), which are known in 
closed form, are required for the calculation of I+(p, q) and i_ (p, q). 

In the case of I±(p, q) where q = p + 2n — 1 the sum S-(2p + 2n — 1) is 
required but this in not generally known in closed form. When q = p + 2n 
where S-(2p + 2n) is required, the integral I±(p,2p + 2n) is calculable in closed 
form. 



6.1 Connection with Nielsen's polylogarithm functions 

The alternating Euler sums S-(r) and the Jordan sums J\(r) and Ji(r) can be 
expressed in terms of Nielsen's generalized polylogarithm functions S n ^ p (z) 

(-lf+f- 1 f 1 ln"- 1 ^) \n p (l-zx) 



S„, p {z) = — — j- / dx, 

(n - iy.pl J x 

and the three special cases £>„ ! (z) — Li n+ i(z), <5>„ jP (1) = s rlj p and S UtP (—1) 
a np where the latter two integrals are given by 

_ (_i)»+p-i /- 1 ln"' 1 ^) ln p (l -x) A 
(n - l)!p! Jo a; 

and 

(-1)"+?-! /- 1 ln"- 1 ^) \n p (l + x) , 
ax. 



(n — l)!p! Jo a; 

In the case of S_ (r) Coffey [13] has given the following integral representation 
forSL(r) 

5 » = TTT /' f" - Tl - ) ln ( X + x ) dx ' 

"(r) J l + ar/ 

which can be rewritten in terms of the polylogarithm function and Nielsen's 
generalized polylogarithm function <S„ : p (z) i.e. 

Sl(r) = Li r+l (-1) +5 r _i, 2 (-1) = (2~ r - l)C(r + 1) + 5v_i, 2 . 

Both Lewin and Kolbig have commented on their inability to express 
^n.p ( — 1) in closed form for even values of n thus leaving the search for a closed 
form expressions for SL(2n+l) = (2~ 2 ™~ 1 ) £(2n+2)+<72n ,2 as an open question. 

It is interesting that closed form expressions for the sums J7i(r), and Ji(r) 
are not generally known for odd values of r except in the case of the two series 

V [^ + l /2) -f 1/2)] = ^ + £ ln 2 (2) - 1 In*(2) - L k{ l/2), 
2 ^ (2k + 1) 3 5760 24 W 24 V/ 4V 7 h 
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V Mfc + l/2) -W2)] = Z ln(2) c(3) _ J3^_*> ^(2)+ J, l n4(2 ) + ^(1/2), 
(2fc) 3 8 wsvy 5760 24 V y 24 W V ' h 



fe=i 



(12) 

which have been obtained (cf. Appendix Al) from work by Coffey [TJ]. The 
methods used to obtain these closed form expressions do not seem to apply for 
any other odd powers of r for the sums in question. Using the sums above, 
S- (3) is then given in closed form by 



S- (3) = 7 - ln(2) C(3) - ^ 4 - ^ ln 2 (2) + 1 ln 4 (2) + 2 Li 4 (l/2). 

It is also possible to find integral representation (cf. Appendix Al) for 
the sums Ji(r), Jiif) and C{r) which are in turn expressible in terms of the 
quantities s nj m and er nj m . Finally the I integrals written in terms of the Nielsen 
integrals are given by 

p 

(-iy J+(p, g) = Y,(-iy C(M) C(P + 1 + 1 - - C(P + 9 + 1) - «p+9-i, 2, 

/i=2 



(-i) p /±(p, g) = E^ 1 )" (2^ p -"-i) Cb+?+i-M)-(2- p - 9 -i)Cb+?+i)-^ +9 -i, 2, 

p=2 

(-l)f J_(p, g) - 21n(2) (2-r-v - 1) Q(p + q) + 2(1 - 2-^-*) «P + 9 + 1) 

+ J2(~ir(2 1 -^ - i) c(m) (2^ p - ? - l) C(p + g + 1 - m) 

Al=2 

+ (i _ 2-P-9) Sp+q _ X] 2 - C(p + g + 1) - 2M(p + g). 

Written in these terms, the integrals 1+ and I_ are seen to be expressible in 
closed form in all cases whereas closed form expressions for I± are obtainable 
for even p + q. Only in the case where p + q — 3 are closed form expressions 
for I± obtainable by the methods used here. 



7 Logarithmic Integrals 

The i(n,m) integrals The integrals 



i(n,m) = / In™ (a;) ln m (l — x) dx — i(m, n), 
Jo 
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arc related to Nielsen's generalized polylogarithms .s„, p . Integration by parts 
yields the relations 

■ / \ m — 1 \ m — 1 

H n i m i _ (_-\\m+n _ C 1 \M H" ~ 1 ; m ~ Mj _ / i \m+n „ 

rt!m! (n — 1)! (m — u)! ^— ' 

If the quantities i*(n, m) = (— l) m+ " m)/m! n! are introduced into the equa- 
tion above we get 

m— 1 m— 1 

i*(rc, m) = 1 + 2J - 1, m - A*) _ X! Sn ,™-f* • ( 13 ) 

fj,=0 fj,=0 

Writing out the lhs of this equation for the case i*(n, m— 1) and subtracting it 
from Eq. (13) we get Pascal's triangular, three term partial difference equation 

i*(n, m) — i*(n, m — 1) — i*{n — 1, m) = — s n>m . (14) 

This can be solved for i*(n,m) in closed form using the Laplace method of 
generating functions |14j with the initial condition i*(0, m) — 1. The final result 
is 

(-l) m + n i(n,m) = (m+n)!-m-n-(m+n-2)!C(2)-m-n!V ( w ~ I/ + m ~ ^ fliz+i) 

^2 (" - ") ! 

. \ ^ (n — u + m — 1)1 „, . , , (n — ^ + m — u)! 

The integrals s rlj p have been given by Kolbig [K] in closed form. Alterna- 
tively, they can be quickly generated with a computer algebra system such as 
Maple or Mathematica using the Kolbig relation 



Vp = (n-l)!p! 



gn+p-l ( iT(l + a)T(l + 0) 



d^daP 1/3 T(l + « + /?) 



a=0=O 



Using this expression, a number of the i(n, m) integrals have been obtained and 
are displayed in the table below. 

Table 1. i(n, m), 1 < n, m < 3 



i(n, m) 


{-l) m+n \m + n)\ 


7T 2 






C(3) 


C(5) 


^ 2 C(3) 


C(3) 


h,i 


2 


-1/6 














h,2 


-6 


1/3 






2 








h,3 


24 


-1 


-1/15 




-6 








«2,2 


24 


-4/3 


-1/90 




-12 










-120 


6 


1/6 






24 


-2 




«3,3 


720 


-36 


-1 


-23/420 


-216 


-144 


12 


36 
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An alternate and more direct way to compute the i(n, to) integrals is due to a 
suggestion by Milgram [16] i.e. noting the relation 



P(v + l,n+l)= / x»{l-x) u dx, 
Jo 

differentiation n times with respect to fi and m times with respect to v and 
evaluated at /j, = v = give 
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m+7f 



-/3^ + 1,/x + l) 



i(m, n) 



i/— o 



As a final remark, we note that the three term partial difference equation 
Eq. (14) can alternately be viewed as a way to compute the s n _ rn integrals 
given independently computed values of the i(n,m) integrals. 



7.0.1 The h(n,m) integrals 



In a similar way the integrals defined by 

h(n,m) = / ln n (x) ln m (l + x) dx, 
Jo 

can be integrated by parts to give 

n!m! 1 j ^ J (n — 1)! (m — /«)! ^ 

which upon introducing the quantities h*(n, m) = (— l) m +" /j(?t,, m)/n\ to! gives 

m — 1 m—1 



n, m — fi 7 



h*(n, m) = 1 + h*(n — 1, m — /i) + 5>, 

M=0 a»=0 



As shown above, a three term partial difference equation for the quantities 
h*[n, m) can similarly be obtained i.e. 

h*(n, to) — h*(n, m — 1) — h*(n — 1, m) = cr n , m - 

In this case, the method of generating functions together with the initial con- 
ditions h(0, to) = —1 + 2 e m (— ln(2)), where e m (a;) is the truncated exponential 
function, gives the solution 

(-l)"+ m /i(n, to) - -(n + to)! + 2n\ [- ln(2)] m 2 F 2 ([-m, n + 1], [ ], 1/ ln(2)) 

. , (n — f + to — /x)! _ 

+m\n\ > > — M , 

; ; ( n — vy.ym — uy. 
(i—i if—i v ' r ' 
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and where 2F2 is a Gauss generalized hypcrgeometric function which reduces to 
a polynomial for the values of the parameters n, m encountered here. We have 

2 F 2 ([-m, n + 1], [ ], 1/ ln(2)) - - ^Jf—j$ fcj ' 



so that 



^ (m-/x)! /*! 

m n , s. 

1 1 (n — z/ + m — //.)! „ 

+m!n! > > — "f- ov, „. 

^ (n - v)\ (m - u)\ ,M 

[1—1 u — l K ' v r 7 

Values for the first few h(n, in) integrals are given below. 

ft(l,l) = 2-2]n(2)-^, 
h(n,m), n + m = 3 

Ml, 2) = 6-^ + ^-81n(2) + 21n 2 (2), 

M2,l) = -6+y + ^C(3)+41n(2), 
h(n,m) 7 n + m = 4: 

h(l,3) = 24 - 361n(2) + 121n 2 (2) -21n 3 (2) + ln 4 (2)/4-7r 2 /2-7r 4 /15 
+3 C(3)/4 - 7T 2 ln 2 (2)/4 + 21 ln(2) C(3)/4 + 6 Li A {\/2), 

h(2,2) = 24- 241n(2)+41n 2 (2)+ln 4 (2)/3-27r 2 /3-7r 4 /12 
-5 C(3)/2 - 7T 2 ln 2 (2)/3 + 71n(2) C(3) + 8Xi 4 (l/2), 

h(3,l) = 24-121n(2)-71n 4 (2)/120-7r 2 /2-9C(3)/2, 

h(n 7 m) 7 n + m = 5 

h{l,A) = -120+ 1921n(2)-721n 2 (2) + 161n 3 (2)-31n 4 (2) + 27r 2 +47r 4 /15 

-3C(3) +7T 2 ln 2 (2) - 211n(2)C(3) - 24Li 4 (l/2), 
h{2,2>) = -120+ 1441n(2) -36 In 2 (2) +41n 3 (2) - 31n 4 (2)/2 + 41n 5 (2)/5 + 3tt 2 

+23tt 4 /60 + [6 - tt 2 - 631n(2)/2 + 21 ln 2 (2)/2] C(3) - 99C(5)/8 

+(3/2- 21n(2)/3)7r 2 ln 2 (2) + [241n(2) - 36] Li 4 (l/2) + 24Li 5 (l/2), 
h(3,2) = -120 + 961n(2) - 121n 2 (2) - ln 4 (2) + 3tt 2 + 11tt 4 /30 

+ [33/2 - tt 2 - 21 ln(2)] C(3) + 87C(5)/8 + tt 2 ln 2 (2) - 24Li 4 (l/2), 
/i(4,l) = -120 + 481n(2) + 27r 2 + 77T 4 /30+18C(3)+45C(5)/2, 
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7.1 Relations between the s n>p and <t„ iP integrals 

The integrals s n , p and a n%v are linearly related. These relations are worth 
examining since they shed some light on the problem of finding closed form 
expressions for the even values of n in the ct„ jP integrals. To show this we note 
that 



"+™ / r(^)r» \1 _ f 1 ln"(x/l+x))ln m (l/(l +a; )) , f 1 ln m ( a /l+ a ))ln"(l/(l +3: )) , 

-a^ ^ TUI+F) )\ ~ X (i+x) ax + 1 (i+x) ax > 



d" 

where the integral representation for the Beta function 

has been used. Expanding the logarithm terms in the integrals above, together 
with the relations 

1 ln g (x)]n p (l+x) , / a \ f 1 ln q - 1 {x)ln p+1 (l+x) , 
ax = — — tt / aaJ, 



(l+z) 



we find 

n m— 1 

/ -i\n+rn„ \ V i\t/n+m-l-l;\~ . \ " / i sfe (n+m— 1 — fe\~ 

S„ + i im — (, TO _! )CTk+l,n+m-k+ 2^, I n J n+m-fcj 

fc=0 fc=0 

or in its more symmetric form 

n m 
fe=l fc=l 

(15) 

As previously pointed out by Kolbig [17] , we see that for a given value of m + n 
there are (for m + n > 6) more of the a v , M integrals than there are equations, 
given the structure of the relations given above in Eq. (15). Closed form values 
for some of the er„. M integrals are given in the tables below 

3%i = -7r 2 /12 
3% 2 =C(3)/8 
5 2 ,i =-3C(3)/4 



Table 2. cr„ !jU , f + /x = 4 



5 1/, n 




tt 2 In" (2) 


In (2) 


ln(2)C(3) 


Lu{l/2) 


01,3 


-1/90 


-1/24 


1/24 


7/8 


1 


0"2, 2 


-1/48 


-1/12 


1/12 


7/4 


2 


03,1 


-7/720 
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Table 3. ov, M , v + \i = 5 





^C(3) 


C(5) 


InW(3) 


Tr^ln 3 ^) 


ln b (2) 


ln(2)ii 4 (l/2) 


L» 6 (l/2) 


0%4 




1 


-7/16 


1/36 


-1/30 


-1 


-1 


(72,3 


1/12 


33/32 


-7/8 


1/18 


-1/15 


-2 


-2 


(73,2 


1/12 


-29/32 












(74,1 




-15/16 













In the case where f + /i = 6, it is not possible to obtain closed form ex- 
pressions for the corresponding <r u ,n values using these methods. However the 
following relations can be found. 



Cl,5 



2 (72, 4 — (74, 



._1_^6 _ l_ n 2 ln 4 (2) + X_ ln 6 (2) + 7_ ln 3 (2) c(3) 



2 



ln^(2)Lz 4 (l/2) + ln(2)£i 5 (l/2) + Li 6 (l/2), 

^ln 4 (2) + iln 6 (2) + ^ln 3 (2)C(3) 



53.^6 1„2,„4, 



15120 24 

iC 2 (3) + 2 ln 2 (2)^ 4 (l/2) + 41n(2)Li 5 (l/2) + 4Li 6 (l/2), 



2(73,3 - 3(7 4 ,2 = jjnT'' 6 ~ |C (3), 



(75, 1 



31 



30240 



A Appendix 

A.l Jordan sums with odd order r 

It is possible to express the Jordan sum J\(2n + 1) with odd argument as the 
integral 



Ji(2n+1) 



1 



4(2n)! J 



ln 2n (a;) In 



1 + x 
1-x 



1 



1 



1 — X 1 + X 



dx. 



This can be done with the help of an integral representation [TS] for the Psi func- 
tion in the numerator and the usual integral representation for the denominator 
within the summand of Eq. (11a) i.e. 



V>(fc + 1/2) - V>(l/2) 
2(2fc + 1)2™+! 



i f°° r 1 (i - t k 

2W0 ^M-[u + iWulA_ 



t) 



dt. 
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Summation over k followed by integration over t produces the desired integral 
representation for J\(2n + 1). 

In the case where n = 1, contains the integrals 



1 ln 2 (x)ln(l - x) , 
dx 



(i 



1-x 180 
l " J( ' ,lH,l+a;) &=|l n (2)C(3)-^ 7 r 4 



1 1 " \x) m^i -f x) 

j— — - ax = 2 - • 72(1 



jf ln2( f^" X) ^ = + $ ln 2 (2) - I ln 4 (2) - 4Xu(l/2), 

/' ^!MM1±^1 ^ = 4Li 4 (l/2) - £ - ^ ln 2 (2) + i ln 4 (2) + f ln(2) C(3), 
Jo It* 

(16) 

which have been given by Gastmann [21,, Lewin [19] and Coffey [13]. These 
have been used to produce the sums given above in Eq. (12). 

Using the same methods employed above in Eq. (lib), the integral repre- 
sentation for J~ 2 (2n + 1) is given by 



J 2 (2n + l) = ^ J o ln 2 «(x) In (I±|) 



1 — x 1 + x 



dx. 



where we note the difference in signs within the integrands of the J7iand Ji 
integral representations. This integral upon expansion can be rewritten as 



l 7 2 (2n+l) = 3( ^^ 1 ln 2 "(x) 



ln(l + x) lnfl — x) d , . , , , 

~Ft~~ \ 7i In l+x In 1-x 

(1 + x) (1 — X) dx 



dx. 



The first two integrals are related to Nielsen's generalized polylogarithms s n<p 
and a Hi p and we get 



J 2 (2n+ 1) = i[s2n,2 + 0-2n,2] 



4(2n)! 



ln 2n (x) 



dx 



ln(l + x) ln(l - x) 



dx. 



The remaining integration within the J?2 integral representation can be carried 
out using the expansion 



dx 



ln(l + x) ln(l - x) =^[21n(2)-V(A:)+V(fe + l/2)] a; 2 ' " 1 . 



fc=i 



We get 



4(2n)! 



ln 2n (x) 



dx 



ln(l + x) ln(l - x) 



dx = \Yl 



[21n(2)-^(fc) + V(fc + l/2)] 



fe=i 



(2fc) 



2n+l 
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The rhs of this equation can be written in terms of the sums defined above as 

\ [j 2 (2n + 1) - C(2n + 1) + 2- 2n - 2 ({2n + 2)] . 

The combined integrals give Jii^ln + 1) as 

J 2 (2n + 1) = \[s 2 n, 2 + v 2n , 2] - C(2n + 1) + 2- 2 ™- 2 C(2n + 2). 

In a similar way, the sum J7i(2n + 1) can be expressed in terms of Nielsen's 
polylogarithms as follows. Adding the integral representations for J\ and Ji 
we get 

M2n + 1) + J 2 (2n + 1) = ^ J* X ^ [x) m (I±|)_^_. 

Expanding the second logarithmic term results in an integral which is in part 
expressible as s 2rii 2 i-e. 

fl 1 2n 



. . , , r ln z "(a;)lni 
^ 1 (2n+l) + J 2 (2n + l) = i S2 „, 2 +2 ( i iyT y o ^— 

Expansion of ln(l + x) / (1 — x) in a power series gives 

i^-±^ = f] [ln(2) + (-l)^ 1 Mk/2 + 1) - <P(h + 1) + ln(2)}] x fe . 
[i X) k=i 

The remaining integral is then given by 

1 f 1 ln 2 "(a;)ln(l + a;) _ 1 [ln(2) + {^(fc/2 + 1) - jj(k + 1) + ln(2)}] 

2(2n)U l-x X ~ 2 t^ l (fc + l) 2 '^ 1 

As in the case of J2{2n +1) these infinite summations are expressible in terms 
of sums defined above and we get 

Ji(2n+1) = |[s 2 „,2-CT2„,2] -M(2n+1). 
More generally, using the same methods one can show for arbitrary r that 

Ji(r) = 5[s r _i,2 - 5 r _i >2 ] - X(r), 

J 2 (r) = i[ Sr -i, 2 + a,.-i, 2 ]-C(r) + 2-''- 1 C(r + l). 



The sum C(r) can also be written as the integral 

C (r) til f 1 ^-\x)lo g (l-x) 

( r) ~ 2^(r-iy.J Q x(l-x) dX ' 
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using the methods used above in connection with the Jordan sums. This integral 
can then be expressed in terms of the Nielsen integrals as 



C(r) = 2^ft (s r ,l + flr-l, 2) = wtt [C(r + 1) + S r _ 



1,2 



As a result of the latter simplified expression, J%(t) can be written as 

Mr) = |[(l-2- r ) s r _i )2 +^-1,2]- 

In the case of the Milgram sum A4(r), it does not appear to be possible to 
express it in terms of the Nielsen integrals. 

A. 2 Approximate values for the alternating Euler sums 

Since it does not appear to be possible to express the Jordan sums of odd order 
as closed form expressions and thus the corresponding alternating Euler sums, 
we have sought approximations for the later quantities. Beginning with the 
relation 

M^H)-Li p (-i)] r^ P M)m(i-t) 

l-t Jo t V ' 

expansion of Li p (—t)on the rhs of this equation followed by integration with 
respect to t gives the alternating series S- (p) i.e. 



I 



' ' ^iif « = t { -=f m + 1) +7] - s-w. as, 



fc=i 



If within the lhs of Eq. (17) Li p (—t) is expanded in powers of t — 1, we get the 
relation 

r 1 Li p (-t) -Li p (-l) dt= _^l j\ t _ lf -r dt j^ s <J) (2W -p_ 1} c(p _ i}) 

(19) 

with the help of formulas for the derivatives |20j of the polylogarithm functions 
i.e. 

r/ A /./,..: /:■"■ ' 



dt k 



= 5>£> (2 i+ ^ p -i)c(p-j), 

t=l ,- =1 



and where Sl"^ are the Stirling numbers of the first kind [2 2) . Integrating over 
t in Eq. (19) we get 

s-(p) = E £ 5 * (2l+ " p - 1)c(p - 
fc=i ■ 3=1 

If the sum on the right-hand side of this equation is truncated to order k = k t , 
we obtain a value of the alternating Euler sum which can be made accurate to 
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any desired number of decimals places. For example taking p = 5 and k t = 10, 
an approximate value of the alternating sum is 

^(^E^E^ (2- 4 -l)C(5-.), 

fe=l ' j=l 

which when written out explicitly takes on the form 

c (r.\ ~ _ 24387227 _ 358039 2 _ 1968329 4 , 2152309 f-(-y\ , 1874237 l (n\ 
°-y°) ~ 1741824000 11197440 130636800 ^ 3456000 ^ 14515200 ^ ' 

an expression which is accurate to nine decimal places. 

It is hoped that future work in this area will produce exact or more accurate 
values of the Euler alternating sums. 
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